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1  Introduction 


In  1978  the  PRL:UT  acoustic  vertical  line  array  model  was 
developed  under  NORDA  Code  520  sponsorship  to  aid  in  the  analysis  of 
PAR  (programmable  acoustic  recorder)  data.  That  model  computes  the 
response  of  a  vertical  line  array  to  a  stationary  point  cw  source  by 
performing  phase  shift  (time  delay)  beamforming  on  the  element  outputs 
(Foreman,  1980).  Acoustic  fields  at  the  elements  are  computed  using 
normal  mode  theory  rather  than  ray  theory  in  order  to  account  for  the 
effects  of  bottom  interaction  and  the  curvature  and  inhomogeneity  of 
wave  fronts  over  the  array  aperture. 

In  1979  the  Doppler  corrected  vertical  line  array  model  OBEROH 
was  developed  in  order  tn  extend  vertical  line  array  modeling 
capability  to  include  the  effects  of  source  motion.  The  model  is  based 
on  a  normal  mode  theory  of  Doppler  effects  (Hawker,  1979)  for  a 
horizontally  stratified  time  independent  medium  and  a  source  moving 
with  constant  speed  and  direction  along  a  nonradial  path  with  respect 
to  a  fixed  receiver. 

The  computer  programs  which  collectively  form  the  OBEROH  package 
calculate  the  pressure  field  produced  by  a  moving  source  over  an 
extended  period  of  time  for  each  element  in  an  array.  The  receiver 
element  response!  are  summed,  according  to  the  user's  shading  and 
steering  requirements,  to  produce  the  beamformed  array  response.  A 
Fourier  transform  is  then  applied  to  form  a  beam  spectrum.  The  purpose 
of  the  OBEROH  program  package  is  to  compute  and  display  beam  spectra 
of  moving  sources  as  the  source  traverses  one  or  more  tracks. 

The  decision  to  base  OBEROH  on  normal  mode  theory  and  to  include 
source  motion  effects,  rather  than  resorting  to  more  conventional  ray 
theory  and  stationary  source  array  modeling  methods,  was  based  on 
several  considerations.  Normal  mode  theory  readily  includes  sediments 
and  substrate,  while  ray  theory  models  must  generally  rely  on  less 
accurate  plane  wave  reflection  coefficient  representations  of  the 
bottom.  Moreover,  normal  mode  theory  can  account  for  the  differences 
in  Doppler  broadening  which  can  occur  across  the  aperture  of  a 
moderate  length  vertical  line  array.  Ray  theory  calculations  tend  to 
be  inaccurate  near  boundaries,  yet  receiver  arrays  are  often  moored 
near  the  bottom.  Finally,  the  response  of  an  array  over  an  extended 
period  to  a  moving  source  will  be  seen  to  be  significantly  different 
from  the  response  of  an  array  to  a  stationary  source  because  of  range 
averaging  effects.  For  all  of  these  reasons,  OBEROH  is  suited  to  the 
prediction  of  array  behavior  under  experimental  conditions  and  ean  be 
used  to  simulate  array  responses  in  the  absence  of  experimental  data. 
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This  report  summarizes  the  mathematical  description  of  OBERON  and 
describes  the  computer  programs  which  implement  the  model.  Because 
this  report  is  intended  for  general  distribution,  no  attempt  is  made 
here  to  describe  the  use  of  the  computer  programs  which  make  up 
OBERON,  nor  is  a  copy  of  the  computer  code  included. 
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2.0  Ebyilfi&l  fiasil  Of  its.  Model 

Figure  1  shows  the  coordinate  system,  geometry,  and  conventions 
used  in  OBCRON  to  specify  array  positions,  array  configurat ions, 
steering  angles,  and  sources.  The  environment  in  which  OBERON  operates 
is  also  illustrated.  Each  of  the  components  in  this  figure  is 
described  fully  in  the  following  sections. 


2.1  Environment 

The  sources  and  vertical  line  arrays  are  immersed  in  an  ocean 
overlying  fluid  sediment  layers  which  in  turn  overlie  a  (possibly) 
solid  substrate.  All  hydroacoustic  and  geoacoustic  parameters  of  the 
ocean  and  sediments  are  restricted  to  depth  variations  only;  no  range 
dependence  is  permitted.  The  substrate  is  homogeneous. 

The  acoustic  parameters  treated  in  the  model  are  density,  p(z), 
sound  speed,  c(z),  attenuation,  niz)  and,  in  tne  substrate  only, 

shear  speed,  c  ,  and  shear  attenuation,  a  . 

s  s 

Density  is  constant  throughout  the  water  column  and  within  each 
of  the  sediment  layers  and  the  substrate,  although  density 
discontinuities  may  exist  at  layer  interfaces. 

Sound  speeds  are  permitted  to  vary  arbitrarily  as  a  function  of 
depth  in  the  water  column  and  sediment  layers.  Sound  speed 
discontinuities  may  occur  at  layer  interfaces. 

Depth  variation  of  attenuation  in  the  sediment  layers,  with 
discontinuities  at  interfaces,  is  also  allowed,  but  the  water  itself 
is  treated  ae  a  lossless  medium.  This  is  consistent  with  the  low 
frequencies  for  which  the  model  war  designed  and  could  easily  be 
general ized. 

The  substrate  is  character ized  by  its  density,  sound  speed, 
attenuation,  shear  speed,  and  shear  attenuation. 

The  coordinate  system  is  chosen  so  that  the  surface  of  the  ocean 
coincides  with  the  x-y  plane  and  depth  increases  along  the  positive  z 
axis.  The  sea  surface,  the  water-sediment  interface,  the  interfaces 
between  sediment  layers,  and  the  sediment-substrate  interface  are  all 
smooth,  planar,  and  parallel.  In  short,  the  entire  environment  is 
horizontally  stratified  to  ensure  the  applicability  of  ordinary  normal 
mo<Je  theory. 
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VEL 


2*2  Sourcti 


OBERON  computes  the  beam  spectra  result  ins  from  point,  time 
harmonic,  omnidirectional  sources  moving  in  straight  lines  at  constant 
depth  z  with  constant  speed  v. 

S 

The  projection  of  a  source  trajectory  onto  the  surface  of  the 
ocean  defines  the  travel  path.  The  acoustic  field  generated  by  the 
source  is  sampled  at  equal  time  intervals  At  for  a  finite  period  T, 
during  which  time  the  source  travels  a  distance  vT.  The  projection  of 
this  line  segment  on  the  travel  path  defines  the  source  track. 

The  point  or  the*  travel  path  nearest  the  origin  of  the  coordinate 
system  is  termed  the  closest  point  of  approach  (CPA)  of  the  path.  Note 
that  the  source  track,  which  is  e  line  segment,  may  or  may  not  contain 
CPA. 


Ths  aouroe  i »  specified  by  its  angular  frequency  to  ,  speed, 
depth,  track  length  (or  sample  period),  CPA,  and  range  fror  the  center 
of  the  track  to  CPA  (r).  The  track  range  from  CPA  is  a  signed 
quantity;  it  is  negative  if  the  source  is  approaching  CPA  and  positive 
if  the  source  is  receding. 

Time  is  measured  with  respect  to  the  moment  that  the  source,  at  a 
distance  r  from  CPA  and  moving  with  speed  v,  would  pass  through  CPA. 
Thus,  time  is  also  a  signed  quantity,  being  negative  for  a  source 
approaching  CPA  and  positive  for  a  receding  source. 


2.3  Vertical  Line  Array  Cunf.l  auras  ion 

The  vertical  line  array  consists  of  one  or  more  point, 
omnidirectional  receiver  elements  distributed  along  the  positive  z 
axis.  The  array  is  specified  by  the  depths  of  the  receiver  elements. 
Spacing  between  the  elements  is  arbitrary  and  elements  may  be  located 
within  the  sediment. 
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3.0  Mathematical  Description 


The  final  product  of  the  beamforming  programs  of  OBERON  is  a 
power  spectrum  for  each  source  track  and  steering  angle.  The 
computations  leading  to  these  results  can  be  divided  into  several 
stages:  (1)  computation  of  the  time  dependent  acoustic  field  at  each 
array  clement;  (2)  calculation  of  the  array  response  (beamforming)  at 
each  time  sample;  and  (3)  production  of  the  power  spectrum  by 
application  of  a  complex  fast  Fourier  transform  to  the  array  response. 


3.1  Calculation  of  Doppler  Corrected  lime  Ssr.ics 

The  derivation  of  the  expression  for  the  Doppler  corrected  time 
series  is  described  in  detail  by  Haujker  (1979).  The  results  are 
summarized  here. 


by 


The  time  dependent  acoustic  field  in  a  lossless  medium  is  given 


Sntr>t) 


U  (Z  )  U  ( Z ) 

nv  s'  nv 


(knR) 
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(r, 


t) 
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£ (0)  =  1 - sin  6 
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where  v^G  is  the  group  velocity  of  normal  mode  n,  un(z)  is  the 
eigenfunction  corresponding  to  the  eigenvalue  kn,  R  is  the  range  from 
the  source  position  on  the  track  to  the  coordinate  system  origin  at 
time  t,  ard  e  is  the  track  bearing  angle  (see  Fig.  1).  The  principal 
assumptions  made  in  the  derivation  are  that  (1)  the  environment  is 
horizontally  stratified,  (2)  the  source  is  moving  in  a  straight  line 
at  constant  speed,  (3)  k  R  >>  1,  and  (4)  v/  vnG  <<  1. 
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The  rang*  invariance  constraint  it  imposed  because  the  Doppler 
corrected  time  series  computation  is  based  on  ordinary  normal  mode 
theory. 

Assuming  that  a(z)  <<  <*>/c(z),  the  loss  due  to  absorption  is 

accounted  for  to  first  order  by  adding  a  small  perturbation  term  i  $ 
to  the  eigenvalue,  where  n 
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where  is  the  depth  of  layer  i  (see  Fig.  1),  Yg  and  Yn  are  first 
order  shear  and  compress ional  perturbation  factors,  and  '*  is  the 
number  of  sediment  layers.  Retaining  only  first  order  perturbation 
terms,  Eq.  (1)  can  be  written  to  include  absorption: 


nz.o 


ip0(2") 


1/2  -  i  (o>t  ♦n/4) 


Vr>t) 


-6  R 
n 


(9) 


Repeated  evaluation  of  Eq.  (S)  at  array  element  depth  z  at  fixed 
time  intervals  as  the  source  moves  along  the  track  yields  the  complex 
time  series  for  that  element. 

The  normal  mode  quantities  un(z)  end  k  n  (Eqs.  (2),  (6),  and 
(9))  are  obtained  from  the  solutions  to  the  depth  separated  wave 
equation 


d  Un(z)  2  2 

- V-  +  (k  (z)-k:)  u  (z)  -  0 

dz  n  n 


(10) 
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A  description  of  the  normal  mode  solution,  and  the  calculation  of 
vG.,  y  •#  and  y  la  given  by  Gonzalez  (I960). 

n  8  n 


3.2  Definition  q£  fiCEjy  Response 

The  response  of  the  vertical  line  array  at  tlise  t  is  the  susi  of 
each  of  the  array  element  responses  at  timr  t  weighted  by  their 
respect lve  shading  and  steering  coefficients: 


R(*,t)  =  E 


m 


s 

m 


w  (♦)  ^(2 

m  n 


(ID 


where  4»  Is  the  steering  elevation/depression  angle,  an  are  real 
amplitude  shading  coaff iciants,  and  the  w  are  complex  steering 
coefficients.  “ 

The  amplitude  shading  coefficients  sm  are  assumed  to  be  equal 
unless  the  user  supplies  other  values  in  an  attempt  to  reduce 
sidelobes.  In  any  event,  the  coefficients  are  rescaled  if  necessary  so 
that  their  sum  is  unity. 

The  plane  wave  steering  coefficients  v  are  given  by  Urick 
(1967),  m 


where 


-ik(z  -z  )sin  * 
m  o 

w  -  e 
m 

k  =  2wf/c 


(12) 

(13) 


3.3  Calculation  of  Beam  Spectra 

If  the  array  response  is  calculated  as  shown  in  the  preceding 
section  for  each  of  N  time  samples  at  constant  time  intervals  At,  the 
result  is  an  array  response  time  series.  The  array  response  time 
series  thus  formed  may  then  be  transformed  using  a  standard  complex 
fast  Fourier  transform  (FFT)  routine  (Brigham,  1974)..  The  bandwidth  of 
the  resulting  spectrum  is  1/At  and  the  frequency  resolution  is 
1/NAt. 


The  side  lobes  appear  if  ;  <n  the  array  response  spectrum  may  be 
•impressed  in  OEERON  by  the  optional  application  of  a  Hanning  complex 
digital  filter  (Brigham,  1974)  at  the  cost  of  reduced  frequency 
resolution. 
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4.0  Software  BeagnAfliinp 


OBERON  Is  a  collection  of  computer  programs  that  perform  the 
three  step  process,  described  in  Sec.  3.0,  under  control  of  the  ARLzUT 
model  operating  system  (Payne,  1900).  These  programs  communicate 
through  disk  files.  There  are  five  main  programs  in  the  OBERON 
package;  MT-ESIS,  RANGDOP,  DOPFMT,  BSPEC,  and  DOPPLT.  The  roles  of 
these  programs  are  discussed  in  the  following  sections.  The 
organization  of  the  program  package  is  illustrated  in  Fig.  2. 

The  modular  structure  of  OBERON  promotes  flexibility  of 
application  for  the  user  and  simplifies  modification  of  the  model.  For 
example  the  time  series  generator,  RANGDOP,  could  be  replaced  with 
another  program  which  computes  the  field  differently,  or  a  measured  or 
computed  ambient  noise  field  could  be  added  to  the  cw  source  generated 
field  without  altering  the  OBERON  package. 


4.1  Time  Series  Generation 

Before  a  time  series  can  be  generated  the  ARL:UT  normal  mode 
modal  NEMESIS  is  run  with  a  user  supplied  hydroacoustic  and 
geoacoustic  description  of  the  environment  and  the  source  frequency. 

It  produces  a  computer  disk  file  containing  3<n,  uR(z),  v®'.,  Yg, 
and  Yn  for  each  normal  mode  n.  Gonzalez  (1900)  describes  this  model. 

RANGDOP  is  then  run  to  compute  the  complex  time  dependent 
acoustic  field  f(z,t).  Figure  3  is  a  flow  diagram  of  the  procedure. 
The  program  reads  the  disk  file  produced  by  NEf'ESIS  and  the  user 
supplied  values  of  attenuation,  source  speed,  source  depth,  range  to 
CPA,  array  element  depths,  number  of  points  to  calculate  in  the  time 
series  (M),  the  rangi*  of  the  center  of  the  track  from  CPA  (r),  and  the 
sample  rate.  Using  Eg.  (9)  the  field  is  computed  from  time  tc  -  M*At/Z 
to  tc  +  MAt/2  at  At  intervals  for  each  source,  array  element, and 
track,  where  tc  ■  r/v.  The  results  are  written  on  a  disk  file. 

Computing  the  time  series  can  be  expensive  because  of  the  large 
number  of  calculations  involved.  Three  measures  were  taken  to  reduce 
the  number  of  calculations.  The  first  of  these  was  the  substitution  of 
•  very  rapid  table  lookup  and  linear  interpolation  procedure  for  the 
standard  FORTRAN  sine  and  cosine  routines. 

The  second  was  the  inclusion  in  RANGDOP  of  a  discrimination 
process  for  identifying  terms  in  Eq.  (1)  which  would  make 
insignificant  contributions  to  the  field  and  deleting  such  terms  from 
the  summation  of  Eq.  (1).  A  term  S  (r,t)  can  be  neglected  if  and  only 
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FIGURE  2  AS-8&81 

SCHEMATIC  ILLUSTRATION  OF  THE  ORGANIZATION  OF  THE  MODEL  OBERON  2  ■  s’ m 
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FIGURE  3 

FLOW  DIAGRAM  OF  THE  TIME  iERIES  GENERATOR  RANGDCP 
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if  its  magnitude  does  not  exceed  a  small  positive  fraction  e  of  the 
magnitude  of  the  term  having  the  greatest  magnitude  for  a  given  range 
R|  that  is 

|Sn(r,t:3 1  *  \Sm(r,t)\K  ,  (14) 
where  |Sm(r,t) |  -max  |Sn(r,t) |  .  (15) 


Substituting  the  expression  ter  Sn(r,t)iEq.  2),  the  test  for  term 
omits  ion  is  written 


-R6 
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A  e 
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<  A  e 
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where 


A  = 
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u  (z  )  u  ( z ) 
n  s  n 

1/2 


(17) 


This  form  of  the  test  requires  that  an  exponential  be  calculated  at 
each  range.  By  taking  logarithms  the  test  is  cast  in  a  form  which 
requires  only  arithmetic  calculations  at  each  range: 


a  -  R5  <  t(R)  +  £ne  >  (1®) 

n  n 

r(R)  =  a  -  R6  »  iU“  (1£0 

v  1  m  n 

a  =  £n  A  .  (2°) 

n  n 


The  last  efficiency  improvement  was  the  development  of  a  time 
series  reformatting  program  DOPTMT.  This  program  reformats  the  time 
series  produced  by  RANGDOP  into  a  form  accessible  to  the  beamformers. 
With  this  program  a  time  series  several  hours  long  can  be  split  into 
separate  source  tracks  of  usei — spec  if ed  duration,  each  with  a 
user-specified  time  delay  between  the  start  of  successive  tracks. 
DOPfTTT  can,  if  requested,  produce  disk  files  of  time  series  generated 
by  overlapping  source  tracks,  yet  no  redundant  time  series 
calculations  are  performed.  Time  series  from  overlapping  source  tracks 
are  often  supplied  as  input  to  the  beamformers.  DOPFMT  can  also  be 
used  to  extract  a  time  series  for  a  track  of  user-specified  length  and 
position  from  the  original,  possibly  very  long  series  generated  by 
RPNGDOP. 
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4.2  Beamformlno 


There  ere  two  beamforming  programs  in  OBERON,  BSPECR  and  BSPECA. 
Both  perform  identical  computations  and  differ  only  in  the  form  of 
output  produced.  The  programs  produce  one  spectrum  for  each  array 
•tearing  angle  and  source  track,  and  create  output  disk  files  suitable 
for  display  or  further  processing. 

The  programs  require  user  input  specifying  the  shading 
coefficients,  steering  angles,  filter  requirements,  and  array  element 
depths.  The  required  time  series  are  read  from  the  disk  file  produced 
by  DOPFTTT. 

BSPECR  produces  spectra  for  one  irray  steering  angle  and  a  series 
of  source  ranges  from  CPA,  that  is,  beam  spectra  at  specified  ranges. 
A  flow  diagram  of  the  procedure  appears  in  Fig.  4.  The  output  can  be 
displayed  using  DOPPLT  as  a  three  dimensional  plot  showing  intensity 
level  (dB)  versus  frequency  and  source  range.  Alternatively,  the  total 
power  in  the  spectrum  can  be  computed  at  each  range  to  form  a  moving 
source  "propagation  loss"  curve  (see  examples). 

BSPECA  computes  the  spectra  at  a  single  source  range  and  a  series 
of  steering  angles.  The  results  can  be  exhibited  as  a  three 
dimensional  plot  showing  intensity  versus  frequency  and  steering 
angle.  Plotting  the  total  power  in  the  spectrum  as  a  function  of 
grazing  angle  yields  a  moving  source  beam  output  curve. 


FIGURE  4 

FLOW  DIAGRAM  OF  THE  BEAMFORMER  BSPECR 
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5  Examples 

Typical  applications  of  the  OBERON  programs  are  illustrated  in 
the  following  examples,  then  appropriate,  results  from  OBERON  are 
compared  and  contrasted  with  results  from  ray  theory  and  with  results 
from  the  stationary  source  vertical  line  array  model. 

The  sound  speed  profile  of  Fig.  5  is  a  typical  deep  ocean 
profile.  The  depth  of  the  water  column  is  4883  m.  The  underlying  200  m 
sediment  is  a  silty  clay  having  a  surficial  sound  speed  of  1523.9 
m/sec,  a  sound  speed  gradient  of  1/sec,  a  density  of  1.42  g/cc  and  an 
attenuation  of  0.1  dB/m/kHz.  The  substrate  has  a  sound  speed  of  5000 
HV’eeo,  a  density  of  2.6  g/cc, and  an  attenuation  of  0.03  dB/m/UHz.  All 
of  the  examples  presented  in  this  section  were  obtained  with  this 
acoustic  environment. 

Normal  modes  and  eigenvalues  were  computed  for  the  acoustic 
environment  of  Fig.  5  and  a  source  frequency  of  25  Hz  using  NEJESIS. 
fill  stationary  and  moving  source  wave  theory  calculations  described 
below  are  based  on  these  results. 

Figure  6  shows  acoustic  field  intensity  versus  depth  in  a  region 
near  the  ocean  bottom  where  a  vertical  line  array  in  deployed.  The 
field  was  computed  assuming  a  stationary  source  100  m  below  the 
surface  at  a  range  of  23.821  km  (the  choice  of  range  will  be  explained 
later).  The  field  is  seen  to  vary  rapidly  over  the  array  aperture. 

The  vertical  line  array  consists  of  seven  uniformly  spaced 
receiver  elements.  The  element  spacing  of  15  m  corresponds  to 
one-quarter  wavelength  at  25  Hz  and  a  nominal  sound  speed  of  1500 
m/sec.  The  array  is  positioned  with  the  deepest  receiver  element  at  a 
depth  of  4710  m,  or  173  m  from  the  bottom. 

The  25  Hz  plane  wave  beam  pattern  for  the  array,  steered 
broadside,  is  displayed  in  Fig.  7.  The  main  lobe  has  a  3  dB  beamwidth 
of  30  deg  and  the  peak  intensities  of  the  side  lobes  are  13  dB  below 
the  peak  intensity  of  the  main  lobe.  The  beamwidth  is  narrow  enough 
and  the  side  lobe  suppression  sufficient  for  the  array  to  distinguish 
between  groups  of  eigenrays  with  different  path  histories  as  the  array 
is  steered  vertically  between  -45  deg  and  45  deg. 

At  this  point  the  array  position  and  configuration  has  been 
established,  the  acoustic  environment  specified,  and  the  normal  modes 
and  eigenvalues  for  that  environment  computed  assuming  a  25  Hz  source. 
Fhere  remains  the  specification  of  the  source  to  satisfy  the  input 
requirements  of  the  time  series  generator,  RANGDOP. 
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The  source  selected  travels  at  10  m/sec  at  a  depth  of  100  m.  Its 
closest  point  of  approach  to  the  array  Is  10  km.  Source  tracks  are 
12 . 775  km  In  length.  The  complex  field  generated  as  the  source 
traverses  a  track  is  sampled  512  times  at  2.5  sec  intervals.  Track 
centers  are  spaced  2.5  km  apart,,  starting  6.4  km  from  CPA  and  ranging 
up  to  106.4  km.  The  tracks  overlap. 

The  next  step  in  the  production  of  beam  spectra  is  the  generation 
of  time  series  for  each  track  and  receiver  element,-  using  RANGDOP  and 
D0Pf?1T. 

Figure  8  shows  output  from  the  beamformer  BSPECR.  Each  line  of 
the  plot  represents  the  array  response  spectrum  for  a  different  source 
track.  The  baseline  of  each  spectrum  represents  a  120  dB  floor  (re  0 
dB  at  1  m)  below  which  the  response  is  ignored.  Baselines  are 
positioned  on  the  plot  so  that  the  range  from  CPU  of  the  track  giving 
rise  to  a  given  spectrum  can  be  read  directly.  Signal  strength  is 
indicated  by  height  above  a  baseline;  the  scale  is  given  by  the  axis 
to  the  right  of  the  plot.  It  should  be  noted  that  the  pressures 
computed  by  RANGDCP  are  scaled  in  the  beamformers  by  iAi  in  order  that 
the  beam  output  levels  not  depend  on  N  (the  number  of  samples)  and  so 
that  the  levels  will  be  consistent  with  stationary  source  propagation 
loss  predictions. 

One  can  follow  the  progress  of  the  source  by  observing  the 
migration  and  evolution  of  the  line  structure  as  the  source  recedes 
from  the  array.  Doppler  shifted  signal  frequencies  decrease 
monoton  leal ly,  approaching  asymptotically  the  limit  of  the  radial 
travel  path.  Signal  strength  decays,  primarily  due  to  geometrical 
spreading  (bottom  reflection  losses  ore  small,  typically  ranging  from 
1  to  4  dB  per  bottom  bounce).  Direct  and  bottom  refracted  paths  show 
the  greatest  signal  strength  and  undergo  the  largest  Doppler  shift; 
paths  which  reflect  from  the  substrate  show  smaller  Doppler  shifts  and 
are  attenuated  while  passing  through  the  sediment.  As  the  source 
recedes,  the  frequency  separation  between  direct  and  reflected  paths 
increases. 

The  beam  spectrum  for  the  track  21.4  km  from  CPA  and  23.621  km 
from  the  array  is  plotted  in  Fig.  9.  Ray  theory  predictions  of  Doppler 
shifts  for  direct,  refracted,  and  reflected  paths  are  indicated.  The 
agreement  is  close  except  when  bottom  image  interference  effects  are 
important. 

A  moving  source  "propagation  loss"  curve  can  be  obtained  from  the 
eeries  of  spectra  shown  in  Fig.  8  by  computing  the  total  energy  in 
each  spectrum  and  plotting  the  result  as  propagation  loss  versus  track 
range  from  the  array.  Such  a  curve  is  compared  in  Fig.  10  with  a 
conventional  propagation  loss  curve  for  a  stationary  source  at  the 
same  series  of  ranges.  Because  the  moving  source  changes  position 
while  the  time  series  is  collected  the  moving  source  propagation  loss 
curve  is  range  averaged  and,  hence,  smoother  than  the  stationary 
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FIGURE  9 

BROADSIDE  BEAM  SPECTRUM  FOR  SOURCE  TRACK  21.4  km  FROM  CPA 
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Figure  11  shows  beam  spectra  as  a  function  of  array  steering 
angle  for  a  track  21.4  km  from  CPA  and  23.621  km  from  the  array.  Uien 
the  total  energy  of  these  spectra  is  plotted  as  moving  source  "beam 
output"  versus  array  steering  angle  another  manifestation  of  range 
averaging  appears.  Figure  12  shows  beam  output  for  the  moving  source 
and  for  three  stationary  sources  positioned  at  the  center  and  both 
ends  of  the  track.  Each  stationary  source  beam  output  curve  contains  a 
null.  The  nulls  occur  at  different  array  steering  angles  for  each 
different  source  position.  No  such  structure  is  discernible  in  the 
moving  source  beam  output  curve  because  the  arrival  structure  varies 
during  the  source  transition  of  the  track. 
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6  Testing 


The  normal  mode  model,  NEJtSIS,  mas  exhaustively  examined  and 
tested  during  its  development  (Gonzalez,  1980). 

The  time  series  generator,  RflNGJOOP,-  has  been  checked  ir.  several 
ways.  The  correctness  of  the  attenuation  calculations  is  indicated  by 
comparisons  oith  an  independently  developed  bottom  loss  model,  as 
discussed  by  Tindls  (1979).  The  tii£«.  series  calculations  wer  tested 
by  computing  the  field  for  a  very  slow  moving  source  (1  ;w.  sec),  as 
well  as  by  checks  of  internal  consistency  and  parametric  dependence  of 
the  output.  Comparisons  of  the  field  at  selected  ranges  with  the  field 
produced  by  a  stationary  source  at  the  sem*  ranges  showed  agreement., 
foalytic  calculations  of  the  Doppler  shift  at  long  ranges  agreed  with 
predictions  produced  by  OBERON.  Ray  theory  predictions  of  Doppler 
broadening  also  agreed  with  those  of  OBERON,  within  the  limitations  of 
ray  theory. 

Using  the  1  mm/sec  source  calculations,  the  output  from  the 
beamforming  programs  was  compared  with  the  output  from  the  stationary 
source  vertical  line  array  model  VLftM  (Foreman,  1980).  Computations  of 
power  at  selected  ranges  showed  agreement  with  a  normal  mode 
propagation  loss  model. 
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7  Conclusions 


A  moving  source  vertical  line  array  propagation  model  (OBERON) 
based  on  a  normal  mode  computation  of  the  acoustic  field  has  been 
designed,  programmed,  tested,  and  made  available  for  use.  The  model 
provides  the  capability  for  computing  vertical  line  array  responses 
which  Include  the  effects  of  source  motion  and  layer  boundaries. 
Possibilities  for  future  extensions  include  adding  measured  (synoptic 
or  archival)  or  modeled  ambient  noise  to  obtain  a  hybrid  "model"  of 
array  response. 
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